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Abstract 

We present the results of a recent novel application of the 'hard exclusive' perturbative formalism to the process 

gg — » MM at large angles, where M (M) is a light meson (anti-meson). As well as discussing the important the- 
Qoretical features of the relevant leading-order gg — > qq{gg)qq{gg) 6-parton amplitudes, we also comment on their phe- 
<*<* nomenological implications. In particular, we consider the central exclusive production of meson pairs at comparatively 

large transverse momentum k±, which is expected to proceed via this mechanism. We show that this leads to various 
^* non-trivial predictions for a range of exclusive processes, and that the cross sections for the 77' and r\ mesons display 

significant sensitivity to any valence gg component of the meson wavefunctions. 



^p-i'l. Introduction 

qj , Central exclusive production (CEP) processes of the 
^ type 

PP(P) -^p + X+p(p) , (1) 
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can significantly extend the physics programme at hadron 
colliders. Here X represents a system of invariant mass 
Mx , and the '+' signs denote the presence of large rapid- 
ity gaps. Such reactions represent an experimentally very 
clean signal and provide a very promising way to investi- 
gate both QCD dynamics and new physics in hadron colli- 
sions, and consequently they have been widely discussed in 
the literature: we refer the reader to [l|, [2|, y, |4J for reviews 
and further references. 

One important feature of such reactions, which we will 
make use of later on, is the dynamical 'Jf = ++ selec- 
tion rule': that is, the produced state must have even C 
parity, while the production of states with non-jj 3 = + 
quantum numbers (where J z is the projection of the ob- 
ject angular momentum on the beam axis) is strongly sup- 
pressed [a, |6j ■ This suppression is predicted in the Durham 
pQCD-based model to be at the level of two orders of mag- 
nitude. 

A particularly interesting example, which has been the 
topic of recent investigations in [3, la|, is the production 
of light meson pairs (X = nn, KK, pp, rj(')r)(')...) at suffi- 
ciently high transverse momentum k± that a perturbative 
approach, applying the Durham pQCD-based model of 
CEP (see e.g. [J] and references therein) and the 'hard ex- 
clusive' formalism (see e.g. [9|, llO[), to evaluate the meson 
production subprocess, can be taken. 

There are two principle, related, reasons to look at such 
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Figure 1: The perturbative mechanism for the exclusive process pp ■ 
p + X + p. 



reactions. First, the parton-level helicity amplitudes rele- 
vant to the CEP of meson pairs, gg — > MM (where M, M 
is a meson, anti-meson) exhibit some remarkable theoreti- 
cal features. Such 'exclusive' gg — > qq{gg)qq{gg) 6-parton 
amplitudes with fixed helicities of the incoming gluons are 
not relevant in a typical high-multiplicity inclusive pro- 
cess, and consequently they have not been studied in this 
context before. Second, they are of phenomenological in- 
terest, being experimentally realistic observables at hadron 
colliders from the Tevatron and RHIC to the LHC, with 
the wide variety of meson states available offering various 
channels with which to probe the non-trivial theory pre- 
dictions for these different processes. We will discuss these 
two aspects in turn. 
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Figure 2: Representative Feynman diagram for the gg — > qqqq(gg) 
process. 
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Figure 3: Representative 'ladder type' diagram, which contributes 
for flavour-singlet mesons. 



2. Properties of the gg — >■ MM amplitudes 

The leading order amplitudes are calculated using an 
extension of the 'hard exclusive' formalism, see [9l Il0| . 
The basic idea is that the hadron-level amplitude can be 
written as a convolution of a (perturbatively calculable) 
parton-level amplitude, T, and a 'distribution amplitude' 
4>, which contains all the (non-perturbative) information 
about the binding of the partons in the meson. The gg — > 
MM amplitude can be written as 
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Axdy(j) M {x)(j>M{y)Txx'{x,y- 1 s,e) , (2) 



where vi is the MM invariant mass, x, y are the meson 
momentum fractions carried by the partons and 6 is the 
scattering angle in the gg cms frame. T\\> is the hard scat- 
tering amplitude for the parton-level process 
gg — > qq(gg) qq{gg), where each qq or gg pair is collinear 
with the meson momentuuLj and has the appropriate colour, 
spin, and flavour content projected out to form the parent 
meson. A, A' are the gluon helicities: for our consider- 
ations there are two independent helicity configurations, 
(±±) and (±t)j which correspond to the incoming gluons 
being in a J z = and \J Z \ = 2 state, respectively, along 
the incoming gg direction. A representative diagram for 
purely qq valence components is shown in Fig. [2] Provided 
the meson k±_ is large enough, all intermediate quark and 
gluon propagators will be far off-shell and the amplitude 
can be calculated using the standard tools of pQCD. 

Originally, in [3, [lQj , the simpler case of initial-state 
photons, 77 — > MM, was considered. In [7|,|8| we extended 
these results to the two-gluon case. Taking the simplest 
case of scalar flavour-non-singlet (tttt, KK...) production, 
this proceeds via the type of diagram shown in Fig. [5] We 
will present the results first, before commenting on them 
at the end of the section. There are 31 Feynman diagrams 
which contribute to the leading-order amplitude, and after 



x For a meson produced with large momentum, \k\, we can to 
good approximation neglect the transverse component of the parton 
momentum, q, with respect to k. 



an explicit calculation we find 
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where 'gg' indicates that the final-state partons are qq 
pairs, '^4,-B' are the gluon colour indices and 



a= (l-x)(l -y)+xy , 
b= (1 -*)(1 -y)-xy . 



(5) 
(6) 



As well as the configuration shown in Fig. [51 the outgoing 
qq pairs can also combine in a second way, with the qq 
pair forming each meson being connected by a quark line. 
A representative such 'ladder-type' diagram is shown in 
Fig. [31 as each meson couples individually to two (isosin- 
glet) gluons, only flavour-singlet states can be produced 
via such a diagram (for e.g. a 7r ± state, this is clear, while 
for a 7T° the uu and dd components of the flavour Fock state 
interfere destructively), while for the case of initial-state 
photons it is forbidden in all cases by colour conservation. 
There are 8 Feynman diagrams which contribute to the 
amplitude, and we find that these give 
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where the label l S,qq' 1 is used to distinguish these ampli- 
tudes, which only contribute for flavour-singlet mesons, 
from ([3l [4]) , which contribute for both flavour-singlet and 
non-singlet states. 

As well as having valence qq components, it is well 
known that the dominantly flavour-singlet rj' (and also, 
through mixing, rf) mesons should have a valence gg com- 
ponent, which also carries flavour-singlet quantum num- 



bers. While different determinations of the rj-rj' mixing 
parameters are generally consistent, the size of the glu- 
onic content of the rf (rj) mesons has remained uncertain 
for many years, with different fits giving conflicting re- 
sults, see for example [ill, [l2( (as well as [8| and references 
ther ein) for a discussion of the theoretical uncertainties 
and [13| for a review of the experimental situation. With 
this in mind, we can extend the above calculation to in- 
clude the case that one or both outgoing qq pairs forming 
the meson states are replaced by gg pairs in a pseudoscalar 
state. Such diagrams will contribute to the gg —> T)(')t)(') 
processes in the presence of any non-zero gg valence com- 
ponent. We find 
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T s + % = T a J_ - 4 -j^t-^x - l){2y - 1) • Tff , (10) 

with 76 Feynman diagrams contributing to (JSJ and 130 
to (ITU)) . We can find no simple form in the case of the 
I J* I = 2 (±,=f) amplitudes, but the angular dependences 
are found numerically to be similar. 

We can see that these amplitudes have some remark- 
able features, which we summarize below: 

1. The quark amplitude corresponding to the case that 
the qq pair forming each meson are not connected 
by a quark line, as in Fig. [21 completely vanishes for 
the case that the incoming gluons are in a J z = 
state, see ((3]). This non-trivial result follows from 
an overall cancellation between the many different 
non-zero contributing Feynman amplitudes, and is 
a direct generalization of the case of the equivalent 
77 — > MM amplitudes, which vanish for neutral 
mesons, when the photons are in a J z = state [9(. 

2. The quark amplitude ([7]) corresponding to the case 
that the qq pair forming each meson is connected 
by a quark line, as in Fig. [31 and which only con- 
tributes for flavour-singlet mesons, does not vanish 
for the case of J z — incoming gluons. Thus we 
observe, from this simple fact, a strikingly different 
behaviour between the cases of flavour singlet and 
non-singlet meson pair production. As we will see, 
this has important phenomenological consequences 
for meson pair CEP. 

3. If we replace the qq pairs forming the flavour-singlet 
meson states by gluons, then again the amplitudes 
([9j [TOj) for J z — incoming gluons do not vanish. 
Moreover, they are in fact identical to the ampli- 
tude ([7]) for a purely valence quark component in 
the mesons, up to overall colour and normalization 
factors. That is, they are predicted to have the same 
angular dependence in the incoming gg rest-frame. 
We emphasize that in these cases, all diagrams al- 
lowed by colour conservation contribute to the total 
amplitude, and not just diagrams of this ladder type 



shown in Fig. [3j for each process ( l S, qq\ 'gg' and 
'gg') the final amplitude receives contributions from 
a set of diagrams which are completely distinct and 
apparently unrelated between the three cases. That 
these final results should be so remarkably similar in 
form is therefore very surprising. 

4. The amplitude ([4]), corresponding to the diagram of 
the type shown in Fig. |3J for | J a | = 2 incoming glu- 
ons, vanishes for a particular value of cos 2 6. This 
behaviour, which at first sight may appear quite un- 
usual, is in fact not completely unexpected: the van- 
ishing of a Born amplitude for the radiation of mass- 
less gauge bosons, for a certain configuration of the 
final state particles, is a known effect, usually labeled 
a 'radiation zero', see for instance [14J,|i5|- The posi- 
tion of the zero is determined by an interplay of both 
the internal (in the present case, colour) and space- 
time (the particle 4-momenta) variables, as can be 
seen in (Q|, where the position of the zero depends 
on the choice of meson wavefunction, 4>{x), through 
the variables a and b, as well as on the QCD colour 
factors. 

5. The parton-level amplitudes corresponding to the 
gg — > MM process can also be understood by a 
novel and instructive application of the MHV for- 
malism (see for example [16[ for a review and [17J for 
a more recent reference). In particular, these ampli- 
tudes for J z — gluons are MHV, so by taking the 
known results for the general 6-parton MHV ampli- 
tudes and substituting in the kinematic and colour 
information relevant to the specific two-meson pro- 
duction process, this allows the previous results to be 
derived in a relatively simple way, without resorting 
to a Feynman diagram calculation. In this way, the 
vanishing of the valence quark amplitude ([3]) for the 
case of J z = incoming gluons follows from a few 
lines of calculation. The identical form of the three 
amplitudes (0 [9j [10]) may also be explained in the 
MHV framework by the fact that the same external 
parton orderings in the MHV amplitudes contribute 
in all three cases. The \J Z \ = 2 amplitude (QJ, for 
example, may also be derived using the 'MHV rules', 
see [l7j for references. We refer the reader to [7J, l8| 
for more technical details. 

6. Although we do not show them explicitly here, the 
valence-quark helicity amplitudes for vector mesons 
(pp, ujuj, (f>4>) have also been calculated in [7]. As in 
(J2J) , for the case that the qq pair forming each meson 
are not connected by a quark line, the J z = ampli- 
tude vanishes. Transversely polarized A = ±1 vector 
mesons cannot be produced within this formalism 
via the diagram of the type shown in Fig. [31 due to 
helicity conservation along the quark line, while an 
explicit calculation shows that the amplitude van- 
ishes for the case of longitudinal polarizations. 



3. Phenomenological implications 

We recall that the CEP of meson pairs at sufficiently 
high k± can be modelled using the Durham pQCD-based 
model for the exclusive production of a generic object X, 
produced via the gg — > X subprocess, see e.g. [l|, la, |6|. 
The case we are interested in is X = MM, for which 
the amplitudes calculated in the previous section can be 
used. We list below the most important phenomenologi- 
cal implications of this in the context of meson pair CEP, 
before briefly discussing the relevance of a separate 'non- 
perturbative' picture in the lower k± region, and the im- 
plications this has for the existing exclusive data. 

1 . We have seen from © that the quark amplitudes cor- 
responding to the case that the qq pair forming each 
meson are not connected by a quark line, as in Fig. [5J 
completely vanishes for the case that the incoming 
gluons are in a J z = state. At leading order, the 
production of flavour non-singlet meson pairs, which 
must proceed via these diagrams, can therefore only 
occur when the incoming gluons are in a \J Z \ = 2 
state. Recalling the L J Z = 0' selection rule discussed 
in the introduction, which strongly disfavours such a 
configuration, this gives the highly non-trivial pre- 
diction that the CEP of flavour non-singlet meson 
pairs (7T7T, KK '...) in the perturbative regime, will be 
strongly suppressed^- We note that e.g. the it + tt~ 
cross section is simply predicted to be twice the neu- 
tral 7r°7r° cross section, due to the non-identity of 
the final state particles. 

2. The above result has an important implication for 
the case of 77 CEP, proceeding via the quark box 
gg — > 77 diagram, see for example [18J]. In general 
the background from exclusive 7r°7r° — > 47 produc- 
tion, when one photon from each 7r° decay is unde- 
tected or the two photons merge, may be important: 
we would naively expect the cross section for such a 
purely QCD process to be larger than in the direct 77 
case. However, if we consider for example the cross 
section at y/s = 1.96 TeV, for the pions/direct pho- 
tons restricted to have E± > 2.5 GeV and \rj\ < 1, 
then we predict a(Tr°TT°)/a('-fj) « 1%. Without the 
suppression due to the J z = vanishing of ©, we 
would instead expect the cross sections to be roughly 
equal. In [19j the observation of 43 77 events at 
yfs = 1.96 TeV with |r?(7)| < 1.0 and E T {l) > 2.5 
GeV, with no other particles detected in —7.4 < 
77 < 7.4 was reported by CDF, and the contami- 
nation caused by 7r°7r° production was determined 
experimentally to be very small, and consistent with 
zero (finding V(7r°7r°)/V(77) < 0.35, at 95% C.L.), 



which provides good support for our approach. An 
observation of 7r°7r° CEP which may hopefully come 
with the increased statistics that a further analysis 
of the existing CDF data can bring [20|, |21| , would 
certainly represent an interesting further test of the 
theoretical formalism. 

gffr [pb/GeV], E± > 2.5 GeV, \i] M \ < 1, </s = 1.96 TeV 
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2 The fusing gluons in CEP can to good approximation be treated 
as being on— shell, in which case the gg and beam axis coincide, and 
so the J z of the gluons in the gg rest-frame must be equal to the 
projection of the centrally produced object's angular momentum on 
the beam axis. 



Figure 4: dcr/dMx for meson transverse energy Ej_ > 2.5 GeV 
and pseudorapidity \r]M\ < 1 at y/s = 1.96 TeV for the CEP of 
meson pairs, calculated within the perturbative framework, with no 
gg valence component included for the tj(') mesons. MSTW08LO 
PDFs are used, and the qq distribution amplitude at starting scale 
fio = 1 GeV is given by the CZ form (see [ld|). 



We have seen that the quark amplitudes correspond- 
ing to the case that the qq pair forming each me- 
son are connected by a quark line, as in Fig. [31 do 
not vanish for J z = incoming gluons. We recall 
that these type of diagrams can only contribute for 
flavour-singlet mesons. From this simple observa- 
tion, we then have the highly non-trivial prediction 
that the CEP of flavour-singlet meson pairs, in par- 
ticular r)'r)', will not be suppressed by the J z = 
selection rule, and so it is predicted to be strongly 
enhanced relative to e.g. irir CEP. In the case of ryq' 
and 7777 production, we also expect some enhance- 
ment, although the level of this is dependent on the 
specific rj-rj' mixing parameters that are taken. An 
observation of these T](')r](') CEP cross sections, or of 
the ratio of &(v(')v(')) to a(inr) or 17(77), f° r which 
various uncertainties (due to PDF choice, soft sur- 
vival factors etc) cancel out, would therefore rep- 
resent another interesting test of this formalism: we 
show representative predictions for such ratios in Ta- 
ble [2 with 0,2 (Ho) = corresponding to the case of 
purely qq valence 7j(') components considered so far. 
We emphasize in particular that the rfrf CEP cross 
section is predicted to be very large, and that the 
branching ratios for the p°j and r)ir + ir~ decays of 
the 77' are sizeable, and may be viable channels. In 
Fig. 2] we plot the predicted Mx distribution for a 
range of meson states, and the relative enhancement 



«;fOo) 


-9.5 





9.5 


vWi)/v{wi) 


210 


1300 


1600 


<r{rfrf)l<r{li) 


3.5 


100 


660 


a{ m )/a{^Ti {> ) 


2.7 


12 


66 



Table 1: Ratios of r](')ri{') and it tt , 77 CEP cross sections at 
y/s = 1.96 TeV with MSTW08LO PDFs JH, for a gg distribu- 
tion amplitude with different choices of a^ (/-in = 1 GeV 2 ) , corre- 
sponding to different normalizations of the gg distribution amplitude 
4>g{ x >Q 2 ) <x a 2 (Q 2 )i guided by the fit of |12|. The qq distribution 
amplitude is given by the CZ form (see |10|0. The meson/photons 
are required to have transverse energy Ej_ > 2.5 GeV and pseudo- 
rapidity |r)| < 1. 
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and suppression of the rjr\' and 7171 cross sections, 
respectively, is clear. In principle r)(')r)(') CEP fol- 
lowed by 77C) — > 77 could therefore represent a fur- 
ther background to exclusive 77 production. How- 
ever, an explicit numerical calculation shows that the 
cross sections for rjrj (and also for 7777) CEP after 
branching to the 47 final state are in fact predicted 
to be a small fraction of the direct 77 CEP cross sec- 
tion for the relevant event selection, with a similar 
result holding for the 7777' final state, see [8] for more 
details. 

We recall that as well as having valence qq compo- 
nents, it is well known that the dominantly flavour- 
singlet rf (and also, through mixing, r\) mesons should 
have a valence gg component, which carries flavour- 
singlet quantum numbers. The CEP of rj' and r\ 
meson pairs, in the perturbative regime, represents 
a novel (and complementary) probe of the size of a 
gg component of these mesons. In particular, the 
contribution from a gg valence component of the 77', 
77 mesons enters at the same (leading) order as the 
purely qq contribution in the CEP cross section. An 
explicit numerical calculation shows that any size- 
able gg component of the 77' (and 77) can have a 
strong effect on the CEP cross section, increasing 
(or decreasing) it by up to ~ an order of magnitude, 
depending on the specific size and sign of thegg com- 
ponent. We show this in Fig.[5j taken from [8j, where 
we plot the Mx distribution for X = vj'vj CEP at 
\[s = 1.96 TeV for a band of possible gg components, 
taking the 'Chernyak-Zhitnisky' (CZ) form [23] for 
the quark distribution amplitude in ©. At the LHC, 
we expect the cross section (for the same event se- 
lection) to be roughly a factor of ~ 3-5 larger for 
y/s = 7-14 TeV, with the particle distributions al- 
most unchanged. 

Due to the identical angular dependence of the am- 
plitudes J2J, © and (jTUJ), the effect of including a 
non-zero gg component of the 77' (77) mesons on the 
TjOvC) CEP amplitudes will to first approximation 
be to multiply them by an overall normalization fac- 
tor. The ratio of the different r)(')r}(') cross sections 
are then determined by the mixing parameters (tak- 



Figurc 5: Differential cross section da/dMx for X = 7/7/ production 
at Vi = 1.96 TeV with MSTW08LO PDFs [3, taking the CZ form 
(see [loll 1 for the quark distribution amplitude, and for a band of 
aSp (fig = 1 GeV 2 ) values, corresponding to different normalizations 
of the gg distribution amplitude (/>q(x,Q 2 ) tx a^iQ 2 ), guided by 
the fit of [131 . The mesons are required to have transverse energy 
E± > 2.5 GeV and pseudorapidity |r?| < 1. 
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irrespective of the size of the gg component. In fact, 
the true ratios are expected to deviate somewhat 
from these values, due to the \J Z \ — 2 contribution 
from the amplitude Q, which can in particular be 
important in the case of the 7777 cross section, where 
the suppression of this contribution, due mainly to 
the J z = selection rule, and that of the flavour- 
singlet components, due to the small mixing angle 
0i, can be comparable. 

Finally, we note that, although it has not been dis- 
cussed above, theoretical studies of meson pair CEP in 
fact have a long history, which predate the perturbative 
Durham approach, see 18j and references therein for more 
details an d 1251] for some recent studies. However, as dis- 
cussed in |18l . |26| , at comparatively large meson transverse 
momenta, k±, CEP should be dominated by the pertur- 
bative 2-gluon exchange mechanism discussed above and 
shown in Fig. [1] At lower k± a study of the transition 
region between these 'non-perturbative' (calculated in a 
Regge theory framework, see Fig. 5 of [18[) and 'perturba- 
tive' regimes may be necessary. This was performed in [18( 
for the case of 7171 CEP. For example, for the (E± > 2.5 
GeV, 1 77 1 < 1) meson pair event selection we have consid- 
ered above, the perturbative contribution is expected to 
be dominant in all cases. 

The predicted cross sections are in general much larger 
in the lower k± (Mx) region, for which the perturba- 
tive approach is not necessarily applicable. The first data 
on exclusive 7r + 7r~ production are already available from 



CDF, which has produced a measurement of the dipion 
mass distribution at M n7r < 5.5 GeV at 0.9 and 1.96 
TeV [27J]. In this case a measurement of the pion k±- dis- 
tribution^ corresponding to the same data would be very 
useful, both as a way to probe the (poorly constrained) 
form factor Fm if), for the off-shell intermediate exchange, 
which enters into the non-perturbative models, as well as 
the onset of any perturbative (i.e. power-like) behaviour 
in the higher k± region, which the pQCD-based model dis- 
cussed above predicttQ We also note that we expect 7r + 7r~ 
(and ir+iiT^ij) data to come soon from CMS [28[, with a 
veto applied on any additional particles within their ra- 
pidity coverage, which should contain a sizeable exclusive 
component. 

A further motivation for studying the CEP of meson 
pairs is as a continuum background to exclusive resonant 
Xc — > kit and Xc — > KK production. In the experimentally 
relevant kinematic region (M„.„. ~ M x , k±iir) ~ M x /2), it 
was found in [18| that the non-perturbative contribution 
discussed above is expected to be dominant, but that with 
suitable cuts on the meson k± and rapidity y, it should rep- 
resent a realistic channel to observe Xco production (both 
Xci and Xc2 production are expected to be suppressed by 
the J z — selection rule). Such studies are, for instance, in 
the LHCb programme [29j . We refer the reader to [18|, |30| 
for more details. 



4. Conclusion 

The study of the CEP of meson pairs iMM = 7T7t, KK, 
r)(')r) (')...) has provided a large number of interesting and 
new results, both phenomenological, in terms of their im- 
plications for the CEP cross sections, and theoretical, in 
terms of the many interesting properties of the relevant 
gg — > MM subamplitudes, calculated using a novel appli- 
cation of the 'hard exclusive' formalism. In this note we 
summarize and discuss these results, concentrating on the 
most important theoretical and phenomenological results 
and predictions, without entering into the detailed treat- 
ment that can be found in [7|, |8|. On the experimental 
side, we emphasize that the CEP of meson pairs is an ac- 
tive field, with existing CDF data on tape and future spe- 
cial low-pile-up runs at the LHC offering the possibility of 



3 More specifically, it is the size of the vector difference 
|(kx(vr + ) — kx(7r _ ))|/2 which is important, so that the effect of 
any non— zero k± of the tttt system due to the p± of the outgoing 
protons is subtracted. 

4 A comparison of the k± (and M^tv) distributions at 0.9 and 
1.96 TeV would probe the size of any possible contamination due to 
proton dissociation Q. It would also be interesting to measure the 
distribution in Ay, the rapidity difference between the pions. This 
would give information about the spin (or more generally intercept 
of the Regge trajectory) of the t-channel exchange M* . 

5 While the Xc(o,2) branching ratios to tttt and KK are the same, 
the exclusive continuum background is in fact expected to be some- 
what smaller in the KK case, see [lal . Moreover, this mode may be 
experimentally more favourable, in particular in terms of suppressing 
the non-exclusive (i.e. with proton dissociation) background, due to 
the generally higher particle ID efficiency in the Kaon case. 



observing the processes and testing the many non-trivial 
predictions presented here. 
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